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<^ . The effect of the spontaneous nucleation of black holes in de Sitter space is 

reviewed and the main steps of the calculation in Nucl.Phys.B 582, 313, 2000 
of the one-loop amplitude of this process are summarized. The existence of 
such an effect suggests that de Sitter space is not a ground state of quantum 
gravity with a positive cosmological constant^ 
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The recent observational evidence for a positive cosmological constant has 
renewed interest in the dynamics of de Sitter space. In particular, quantum 
gravity in de Sitter space is attracting new attention (see for example || 
and references therein). In this connection I would like to review in this 
Qh! lecture the effect of semiclassical instability of de Sitter space with respect 

to spontaneous nucleation of black holes. Such an instability suggests that 
de Sitter space is not a ground state of quantum gravity with a positive 
cosmological constant A > 0. 

The existence of the semiclassical instability in quantum gravity at finite 
temperature has been known for a long time. This is a non-perturbative 
effect mediated by certain gravitational instantons. In the A = case the 
effect was described in detail by Gross, Perry and Yaffy [14[], while the case 



with A > was considered by Ginsparg and Perry |[L3|| . The simplest way 
to understand this effect is to consider the semiclassical partition function 
for Euclidean quantum gravity with A > as the sum over gravitational 
instantons. It turns out that in this sum there is one term, determined by 
the S 2 x S 2 solution of the Euclidean Einstein equations, whose contribution 
is purely imaginary, because this (and only this) instanton is not a local 
minimum of the action. This renders the partition function and free energy 
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complex, and the existence of an imaginary part of the free energy signals 
that the system is metastable. Since the S 2 x S 2 instanton is the analytic 
continuation of the extreme Schwarzschild-de Sitter solution, one can argue 
that this instability results in the nucleation of black holes in de Sitter space. 
The energy for this process comes from the de Sitter heat-bath. 

This effect is inherently present in quantum gravity with A > 0, since 
the value of the de Sitter temperature is fixed by A, and this gives rise to 
the thermal instability. For comparison, in the A = case the situation is 
different, since the temperature is then a free parameter and one can take 
the zero temperature limit. In this limit Minkowski space is stable |14 . 



Summarizing, the S 2 x S 2 instanton is analogues to the "bounce" solu- 
tion ||, and so it is responsible for the formation of bubbles of a new phase. 
Its physical negative mode living in the physical TT-sector should not be 
confused with the negative modes from the conformal sector [jl(J||. The latter 



exist for any background, and it seems that they do not imply any physi- 
cal instabilities but emerge merely as a result of a bad choice of variables 
in the path integral. It seems likely that if one was able to start from the 
path integral over the phase space and then make it covariant in some way, 



these modes would be absent [18|. Negative modes from the TT-sector are 
on the other hand physically significant. In fact, having factorized the con- 
formal sector, it is only these remaining modes that render the partition 
function divergent. This reflects the breakdown of the canonical ensemble 
for gravity due to black holes, whose degeneracy factor grows too fast with 
the energy. After complex-rotating the physical negative modes, the partition 
function becomes finite but complex, thus again indicating the breakdown of 
the canonical description and the existence of metastability. 

A very nice feature of the S 2 x S 2 instanton is its high symmetry. This 
allows one to carry out all the calculations exactly (in one loop), which 
is an exceptional situation in quantum gravity for non-trivial backgrounds. 
Such a calculation had not been done though until our recent paper with 
Andreas Wipf where a complete derivation of the amplitude of the 



process was presented. In brief, we succeeded to explicitly determine the 
spectra of all the relevant fluctuation operators and to exactly compute their 
determinants within the (^-function scheme. Since calculations of this type 
are rather involved (and also rare), I would like to briefly summarize below 
the essential steps of our procedure, referring to [l9j for more details and 
references. 



2 



1 Qualitative description 



Let us consider the partition function for the gravitational field 

Z = J Dlg^e- 1 , (1) 
where the integral is taken over Riemannian metrics with the action 

with A > 0. The extrema of this action are solutions (M,gnu) of the Eu- 
clidean Einstein equations 

Rfw = Ag^ . (3) 

Such solutions are called instantons; all of them are compact. In the semi- 
classical approximation, for GA <C 1, the path integral (P is given by the 
sum over all instantons 

'DetA m 



Here 7m is the action of the Z-th instanton and the prefactor comes from the 
Gaussian integration over small fluctuations around this background, A[q 
being the corresponding fluctuation operator. 

The leading contribution to this sum is given by those instantons whose 
action is minimal. One can expect that these will be the instantons whose 
symmetry is maximal. Among all solutions to equations fl3|) there is one 
which is a maximal symmetry space - this is the four-sphere with radius 
3/A and the standard metric which can be written as 

ds 2 = (1 - K 2 ) dt 2 + + r 2 {d$ 2 + sin 2 #V) , (5) 

3 1 — -^r 1 



3 ' 
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with r < y 3/A. The action of this S instanton is I — — ?>it/ KG. Hence, for 
GA <C 1, the path integral ([!]) approximately is 

exp(37r/AG) 

Z ~ Z[S4] - /DetA [g4] - (6) 
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One can rewrite this expression as the partition function of some thermal 
system, Z = e~ F ^ T . Here the inverse temperature 1/T = 2 , kJ^ is determined 
by the proper length of geodesies on S 14 (all of them are periodic), and up to 
subleading terms the free energy is F = ~ 2( /^ - 

One can wonder to which physical system refer these temperature and 
free energy ? In fact, they relate to the de Sitter space ||. Upon analytic 
continuation t —*■ it the metric (|5|) becomes the de Sitter metric restricted 



to the region inside the cosmological horizon, r < r c = y3/A. Let us call 
this region the Hubble region, or the causal diamond. De Sitter space has 
the temperature T determined by the surface gravity of the horizon, the 
entropy S = nr 2 /G determined by the area of the horizon, and the free 
energy F = —TS contained inside the horizon [^, |J. The values of T and 
F exactly coincide with those obtained above from the geometry of the S 4 
instanton. 

Summarizing, the partition function of quantum gravity with A > and 
GA <d 1 is determined by the thermodynamic parameters of de Sitter space, 

Z « Z m = e- F ' T . (7) 

Let us now consider corrections to this formula due to other gravitational 
instantons contributing to (f|): 

i+E^ 1 ). (8) 

For AG <C 1 all terms in this sum are exponentially small and can be ne- 
glected, if only they are real. If there are complex terms, their contribution 
will be physically important despite their smallness. Such complex terms can 
arise due to those instantons which are not local minima but saddle points of 
the action. If the number of their negative modes is odd, the determinants of 
t he corres ponding fluctuation operators will be negative, and the prefactors 
wTJetA^ in @ will therefore be complex. 

The question therefore arises: if among solutions of Eqs. (|3]) there are those 
which are not local minima of the action ? A theorem quoted by Gibbons in 
Ref. 0] states that for A > there exists only one such solution, which is the 
geometrical direct product of S 2 x S* 2 with the standard metric: 

ds 2 = - ( d^ 2 + sin 2 tfx d^ 2 + d$ 2 + sin 2 <& d^j . (9) 
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Its action is / = —2tt/GA. As we will see, the solution has exactly one 
negative mode. In view of this, the partition function can actually be ap- 
proximated by the semiclassical contributions of only two instantons, S" 4 and 
S 2 x S 2 : 

V*s 2 A~ex P { -i (f-T^^]\, (10) 



V Z[si] J [ T \ Z[s*] , / 

where Z[s 2 xs 2 ] is purely imaginary. Rewriting this as Z ta e - ^ 7 the free 
energy T here will have the real part coinciding with the free energy F in 
the Hubble region, and also the exponentially small imaginary part 

The existence of the imaginary part of the free energy indicates that the 
system is metastable JR], [If. The conclusion therefore is that de Sitter space, 



which is classically stable, becomes unstable semiclassically when the non- 
perturbative effects are taken into account. 

The qualitative description of this instability was first given by Ginsparg 
and Perry fl3f| , who argued that the quantum decay of de Sitter space will 
lead to the spontaneous nucleation of black holes. The logic is that, as 
the decay is the tunneling transition mediated by the S 2 x S 2 instanton, 
the structure of the configuration created during such a transition will be 
inherited from that of this instanton. The basic observation is then that 
the S 2 x S 2 instanton can be obtained via the analytic continuation of the 
Lorentzian Schwarzschild-de Sitter solution [0, [13|, [| 

dv^ 

ds 2 = -Ndt 2 + j^ + r 2 (d$ 2 + sin 2 tidy?) . (12) 
Here N = 1 — — — 4^ 2 , and for 9M 2 A < 1 this function has roots at 

r 3 ' 

r = r + > (black hole horizon) and at r = r ++ > r + (cosmological horizon). 
If one analytically continues via t — > it, the metric becomes Euclidean, but 
then one has to restrict the range of r to r + < r < r ++ , since N is negative 
otherwise. In addition, one has to identify t with a suitable period, since 
otherwise the geometry will have conical singularities at r = r + and r = r ++ . 
In general one cannot remove both singularities at the same time, since the 
period of the identification of t is determined by the surface gravity, which 
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is different for the two horizon. However, in the limit r + — > r ++ — > 
the surface gravities will be the same and both conical singularities can be 
removed at the same time. Although one might think that nothing will 
remain of the solution in this limit, this is not so. The limit r + — > r ++ 
implies that 9M 2 A = 1 — 3e 2 with e — > 0. One can introduce new coordinates 
$1 and <f i via cos-#i = (vAr — l)/e + e/6 and = \/Aer. Passing to the 
new coordinates and taking the limit e — > 0, the result is exactly the metric 

The conclusion is that a tunneling transition via the S 2 x S* 2 instanton 
will create an extreme Schwarzschild-de Sitter black hole. The radius of such 
a black hole is equal to the radius of the cosmological horizon, and so it will 
fill completely the Hubble region. However, the total volume of de Sitter 
space is infinite and it contains infinitely many Hubble regions. The black 
holes will emerge in some of these regions, but most of the regions will remain 
empty. The number of the filled regions divided by the number of the empty 
ones is the probability for the black hole nucleation in one region. This is 



proportional to in (IT]). As a result, black holes will nucleate with a 

certain probability all over the space, like the bubbles in boiling water. 

The temperature of the nucleated black holes can be read of from the 
S 2 x S 2 metric as the inverse length of the equator of any of the two spheres: 
Tbh = (the same value can be obtained from the Lorentzian solution 
( 12[) [|J). This is different from the temperature of the de Sitter heat bath. 



The origin of this discrepancy can be traced to the causal structure of de 
Sitter space: the fluctuations cannot absorb energy from and emit energy 
into the whole of de Sitter space, but can only exchange energy with the 
Hubble region. Thus the energy exchange is restricted. As a result, the local 
temperature in the vicinity of a created defect may be different from that of 
the heat bath. 

Using finally the classical formula of Langer and Affleck [TJ|, []J , the qua- 
siclassical decay rate of de Sitter space is 

\Z 2 2 I 

r = 2|S(.T)| = 2T 1 [s2 x52]l , (13) 

Z{s 4 ] 

where T = ^r^f is the temperature of the de Sitter heat bath. This formula 
gives the probability of a black hole nucleation per per Hubble volume and 



per unit time of a freely falling observer [19|. 
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2 Calculation of the path integral 



Let us now briefly describe the one-loop calculation of the partition functions 
Z[ S 2 xS 2] and Z[£4] entering formula (|13|). We use the standard Faddeev- Popov 
procedure and complex-rotate the conformal and TT-negative modes. Only 
the ghost operator has zero modes in its spectrum, which modes are asso- 
ciated with the background isometries. We integrate over these modes non- 
perturbatively, with the result (Eq.([24])) being proportional to the volume of 
the isometry group. The main steps of the procedure are as follows. 

For small fluctuations h^ v around an instanton {M-^g^v) the action ex- 
pands as I[g^ u + h^y] = I[g^ v ] + 5 2 I + . . . . Since 5 2 I has many zero modes 
associated with diffeomorphisms, one has to add a gauge fixing term, which 



is chosen to correspond to the covariant background gauge [|lT 



5 2 i g = 7 (v„h; - 2±1 v p h, vx - 1±± v'/A . (14) 

Here 7 is a real parameter and the scalar product is defined as 

similarly for vectors and scalars. The one-loop partition function then reads 

Z = e- 1 J D[V] 2>fp exp (-5 2 I - 5 2 I g ) , (15) 

where the Faddeev-Popov factor is given by the integral over all diffeomor- 
phisms 

(Pfp)- 1 ^ J D[Q exp (-5 2 I g ). (16) 
It is convenient to use the Hodge decomposition for fluctuations, 

V = (f>,xu + ^ hg^ + V M £„ + V p ^ - - g^V^ , (17) 

where 4> pi/ is the transverse tracefree (TT) part, V M 0^ = 0^ = 0, and h is 
the trace. The longitudinal vector part also decomposes as 

£m = ^ + V„X + C (18) 
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where = and the harmonic piece £^ vanishes for simply-connected 

manifolds. With these decompositions the gauge fixing term and the gauge- 
fixed action S 2 I g f = 5 2 I + 5 2 I g are diagonal: 

5 2 I g = 7<77 M , A?rf > + ±-((h + 2A^ x ),A (h + 2A^)> . (19) 

#hf = \ (<T, + 7<»7m> Aft/*) + ± (x, AoAoA^x) - (/i, A /i> . 
Here the tensor fluctuation operator is 

A 2 </v = -V CT V CT M , - 2R pai/ ^ , (20) 
the vector operator and scalar operators are 

Ai = — VffV 7 - A , A = -V CT V CT , (21) 

while A = 3A - 4A and A^ = 7A0 - A also h = h- 2V jU £". 

To compute the path integrals in (JT^,([nj), all fields are expanded with 
respect to the bases associated with complete sets of eigenstates of the op- 
erators A 2 , Ai and A with some Fourier coefficients C. For example, 
h = J2 n C n h n where A h n = X n h n . The quadratic actions ([0J) then re- 
duces to quadratic forms in the coefficients C, and the path-integration in 
(P~5|) , (p~6|) is performed by integrating over all dC. The perturbative path 
integration measure is defined as the square root of the determinant of the 
metric on the function space of fluctuations: 

D[h, u ] ~ ^Det((dh, u ,dh^)) , D[Q ~ ^Det((d^,d^)) , (22) 

where the differentials refer to the Fourier coefficients C of the expansions of 
h^ u and For example, dh = J2 n dC n h n . The normalization is chosen such 
that fD[h flv ]exp(-& (h^bT)) = 1 and / D[Q exp (-^ p)) = 1, 
where \i Q is an arbitrary renormalization parameter with the dimension of an 
inverse length. 

Using (|TT|),(|i~8|) the metrics on the space of fluctuations are 

(dh^, dlf v ) = d^ v ) + 2(%, Atfrf) + (d X , A A rfx) + - (dh, dh) , 

(d^,de) = (d Vl ^d^) + (dx,A dx}. (23) 
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The path integrals in flI5|),(|IED reduce then to infinite products of ordinary 
integrals over the Fourier coefficients C. Most of these integrals are Gaussian, 
and their computation gives products of eigenvalues a s for coexact vectors 
and eigenvalues e& for TT-tensors in the resulting formula for Z in Eq.([25|) 
below. There are, however, also non-Gaussian integrals, which can be of 
three different types. 

• First, for the S 2 x S 2 instanton there exists a TT tensor mode with 
a negative eigenvalue e_ < 0. The integration over this mode is carried 
out with the complex contour rotation || leading to the complex factor 



^ne g = A t o/(2iy |e_ |) in fl25|). For any other instanton Q neg = 1. 

• Next, for any instanton there are infinitely many conformal negative 
modes associated with the 'wrong' sign term — j|- (h,Aoh) in ([19]). This is 
the manifestation of the well-known problem of conformal sector in Euclidean 



quantum gravity JT0|. Although its real understanding is lacking, it seems 
that this problem is essentially an artifact of the bad choice of variables 
in the path integral [13]. The prescription is then to perform the complex 



rotation h — > ih in the space spanned by positive eigenstates of the operator 
Ao [|BJ. After this the quadratic forms become positive-definite and the 
integrals can be computed. Remarkably, the resulting effect of conformal 
modes is then exactly cancelled by the contribution of the exact parts x °f 
the vectors. There is only one scalar mode which contributes to the final 
answer: the constant conformal mode present for any background and giving 
the factor fl con { = to (p5]). Only for the S A instanton there are additional 
5 conformal Killing scalars which also contribute, and so the answer in the 
S 4 case is fi conf = — fMl x " 



V2A VV 3 Mo, 

• Finally, after fixing the gauge, there remain a finite number of zero 
modes of the vector operator Ai. These are associated with the background 
isometries. Such modes are not contained in the fluctuation measure D]^^], 
since they do not contribute to the metric (dh^^dW) in fl23|) , but they are 
contained in the ghost measure D[^] and should therefore be taken into 
account. In fact, the effect of these modes turns out to be very impor- 
tant. In treating these modes, we follow the approach of t'Hooft |T5[ and 



Osborn |17[], whose idea is to carry out the integration over these modes non- 
perturbatively, which amounts to integrating over the isometry group 7i. The 
corresponding integration measure must be proportional to the Haar mea- 
sure. If TL acts on M. via x M — > x M ((7j), where j = 1, . . . dimH, the vector 
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zero modes are the Killing vectors Kj = -^j- = ffrgf^j and the integration 

over these modes gives frilling = / (jlj ^5= g^T ) dfi(C) . Here is com- 
puted at Cj = and the normalization of the Haar measure dfj,(C) of TC is 
fixed by the condition that at the unity, Cj = 0, the perturbative measure 
is reproduced: dfi(C) — > TTj dCj as Cj — > 0. For the 5 2 x 5 2 instanton the 
isometry group is H = SO (3) x 50(3), while in the 5 4 case TL = 50(5). 
This gives, respectively [pi 



Killing 



647T 4 (/i ) 12 

27(AG) 3 : 



Killing 
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1287T 6 (/i Q ) 



20 



(AG) 5 ' 



(24) 



All this finally leads to the following expression for the one-loop partition 
function of small fluctuations around a given instanton background: 



f2 ne gf2 con f 



Killing 




(25) 



where all dependence on the gauge parameter 7 has disappeared. 

In order to use this formula, we need to explicitly know the spectra of 
the tensor fluctuation operator ( p0[ ) on the space of the TT-tensors subject 
to V M 0^ = 0^ = 0, and also those for the vector operator (|2l| ) on the space 



of coexact vectors V^r^ = 0. For the 5 2 x 5 2 instanton one finds [T^ 



operator 


eigenvalue 


degeneracy 




A 2 


-2A 
2A 

(j(j + l)-2)A 

j(j + 1)A 
(jT(ji + l)+J 2 (j2 + l)-2)A 


1 

9 

2(2j + 1) 
18(2j + 1) 
5(2j 1 + 1)(2j 2 + 1) 


J>2 
J>2 

J1,J2>2 






+ 1)-2)A 
0T(ji + l)+J 2 (j 2 + l)-2)A 


6 

2(2j + 1) 
3(2j 1 + l)(2j 2 + l) 


3i, 32 > 1 


A 


(Ji(ji + 1)+J 2 (j 2 + 1))A 


(2j 1 + l)(2j 2 + l) 


Jl,j2>0 



These spectra contain one negative tensor mode and six vector zero modes 
associated with the 50(3) x 50(3) isometry group; all other modes are 
positive. For the 5 4 instanton one has 
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operator 


eigenvalue 


degeneracy 




A 2 


fj(j + 3) 


|(j-1)(j+4)(2j + 3) 


J>2 


A: 


f(j(j + 3)-4) 


ij(j + 3)(2j + 3) 




A 


fj(j + 3) 


i(j + l)(j + 2)(2j + 3) 


J>0 



In this case all eigenvalues are positive, apart from ten vector zero modes 
which are generators of SO (5). 

The next step is to compute the infinite products of these eigenvalues 
in Eq . (|25|) . This can be done with the use of the (^-function regularization, 
in which scheme the product of an infinite discrete set of numbers a s by 
definition is 

' ' ' ' (26) 



rj- =exp{ _ C '(0)-C(0)ln/x} 

s A 4 



where the (^-function 



(27) 



The spectra in the tables above lead to the ^-functions of the following three 
basic types 



oo oo 



Z(k, v\z) 
C(k,u\z) 
Q(k, u,c\z) 



EE 

n=k 
oo 

E 



(2n+ 1) (2m + 1) 



n=km=k {(2^ + 
(2n + 1 



l) 2 + (2m+l) 2 + z/} 2 



(28) 



n=k 

oo 

E 

j=k 



{(2n+l) 2 + ^ ' 
(2j + 3)(j(j + 3) + c) 
0*0* + 3) + i/}* 



(29) 



where dt(z) must be large enough to ensure convergence. The analytic con- 
tinuation of these expressions to arbitrary values of z has been carried out 
in [IJ|. This gives for z = 



Z(k,i/\0) 

C(k,u\0) 
Q(k,v,c\0) 



32 
1 



v 



24 



i fc V + 2fc 4 --fc 2 + — , 
2 3 360 



-v - k 2 

12 4 

fc 4 



2k 



3 



(30) 

(2c+l)y + (3-2c)fc + ^ + (4-3^-li, 



which quantities define the 'regularized numbers of eigenvalues' of the op- 
erators and determine the anomalous scaling behavior of the determinants. 
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The values of the determinants themselves are given by the derivatives of the 
^-functions at zero, which can be expressed in quadratures and computed 
numerically. The values needed in Eq.((§) are Q'(2,0, -4|0) = 3.72344, 
Q'(2, -4, 0|0) = 6.65246, and Z'(2, — 10|0) = -18.3118 @. 



Collecting everything together, the one-loop partition functions in the 
(-function regularization scheme are 

_98 / Oil \ 

Z [S 2 XS2] = -i 0.3667 x (AG) fi 45 exp (-) (31) 

and 

- -521 / 37T \ 

Z m = 0.0047 x (AG) 5 fio 45 exp ( — ) . (32) 

Up to my knowledge, such closed expressions were obtained for the first time 
|19 . Here the numerical prefactors are determined by the (-regularized 



in 



determinants, the factors of AG come from the background isometries (see 



Eq.(24)). The powers of \i Q are the anomalous dimensions, they receive con- 
tributions from all modes, including the isometries. As is seen from Eq.fl2"4|), 
the anomalous effect of the isometries is considerable, and had it been ne- 
glected the result would be very different. The anomalous dimension — ^ 



for S 4 in fl32|) agrees with the analysis of p|, p0[1 . 
The last step is to insert into Eq. flllf ) to find 



T = 14.338 VX (GA)- 2 (fx A)^ exp (-^) • (33) 

This gives the probability of the spontaneous nucleation of a black hole inside 
the finite volume enclosed by the de Sitter cosmological horizon per unit time 
of a freely falling observer. The formula applies for AG <C 1. Due to non- 
renormalizability of gravity, the cutoff parameter fi Q remains undetermined; 
for estimates one can for example set \i = G. The subsequent real time 
evolution of the created black holes is an issue requiring special study B. 
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